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Abstract 

In this paper, we consider a simple kinetic model of economy involving both ex- 
changes between agents and speculative trading. We show that the kinetic model ad- 
mits non trivial quasi-stationary states with power law tails of Pareto type. In order to 
do this we consider a suitable asymptotic limit of the model yielding a Fokker-Planck 
equation for the distribution of wealth among individuals. For this equation the sta- 
tionary state can be easily derived and shows a Pareto power law tail. Numerical 
results confirm the previous analysis. 

Keywords. Econophysics, Boltzmann equation, wealth and income distributions, 
Pareto distribution. 
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1 Introduction 



Microscopic models of simple market economies have been recently introduced by several 
authors |ni[7l lll[ll4j . The main idea is that an economic system composed by a sufficiently 
large number of agents can be described using the laws of statistical mechanics as it 
happens in a physical system composed of many interacting particles. The details of the 
economical interactions between agents characterize their wealth distribution. 

The study of wealth distributions has a long history going back to more then one 
century ago with the Italian sociologist and economist Vilfredo Pareto which studied the 
distribution of income among people of different western countries and found an inverse 
power law [23|- More precisely if f{w) is the probability density function of agents with 
wealth w we have 



Pareto mistakenly believed that power laws apply to the whole distribution with a universal 
exponent n approximatively equal to 1.5. Later, Mandelbrot, proposed a weak Pareto law 
that applies only to high incomes [17,. 

Income data from developed countries show Pareto-like behaviors with order 1 expo- 
nents for large values of incomes (USA ~ 1.6, Japan ~ 1.8 — 2.2 |llj). It is common, in 
fact, that 90% of the total wealth is owned by only 5% of the population. According to 
data, across the full range of income, we should expect the probability density function 
f{w) to increase at low income, reach a maximum and, finally decrease with increasing 
wealth HI]- 

In our discussion we use the terms distribution of income (money) and distribution 
of wealth interchangeably (although the two distributions may not exactly coincide there 
is a strong dependence on one another). Income, in fact, is only one part of wealth, the 
other part being material wealth. Material products have no conservation law and their 
monetary value is not constant (they can be manufactured, destroyed, consumed, etc.). 
The distinction between the two will not influence the conclusion drawn here from the 
dynamical model of economic interactions we consider. 

Following this line of thought, we consider here a very simple model of an open market 
economy involving both assets exchanges between individuals and speculative trading. In 
our non-stationary economy the total wealth is not conserved due to a random dynamics 
which describes the spontaneous growth or decrease of wealth due to investments in the 
stock market, housing, ect. It is important to note that this mechanism corresponds to 
the effects of an open market economy where the investments cause the total economy 
to growth (more precisely the rich would get richer and the poor would get poorer). 
The conservative exchanges dynamic between individuals redistributes the wealth among 
people. 

We shall show that this kinetic model, introduced in , gives in a suitable asymptotic 
limit (hereafter called continuous trading limit) a partial differential equation of Fokker- 
Planck type for the distribution of money among individuals. 

The equilibrium state of the Fokker-Planck equation can be computed explicitly and 
is of Pareto type, namely it is characterized by a power-law tail for the richest individuals. 
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Our result can be related to the work j^, where a similar equation has been derived by a 
mean-field limit procedure. 

The mathematical methods we use are close to those used in the context of kinetic 
theory for granular flows, where the limit procedure is known as quasi-elastic asymptotics 
[2i5\ . We mention here a similar asymptotic analysis performed on a kinetic model of a 
simple market economy with a constant growth mechanism j25| . 

The paper is organized as follows. In the next section, we introduce the binary inter- 
action between agents, which is at the basis of the kinetic model. The main properties of 
the model are discussed in section 3. These properties justify the continuous trading limit 
procedure, performed in section 4. The limit is illustrated by several numerical results in 
section 5. 

2 A kinetic model of money asset exchanges 

The goal of a kinetic model of simple market economies, is to describe the evolution of the 
distribution of money by means of microscopic interactions among agents or individuals 
which perform exchange of money. Each trade can indeed be interpreted as an interaction 
where a fraction of the money changes hands. We will assume that this wealth after the 
interaction is non negative, which corresponds to impose that no debts are allowed. This 
rule emphasizes the difference between economic interactions, where not all outcomes are 
permitted, and the classical interactions between molecules. 

^From a microscopic view point, the binary interaction is described by the following 
rules (we refer to j22j for more details) 

■w' = {1 — + JW^: + 7JW 

= (1 — 7)u;* + + rj^w^ 

where {w,w^) denote the (positive) money of two arbitrary individuals before the trade 
and {w',w'^) the money after the trade. In (pQ) we will not allow agents to have debts, 
and thus the interaction takes place only if tt;' > and w'^ > 0. In the transaction 
coefficient 7 € (0, 1) is a given constant, while rj and rj^, are random variables with the 
same distribution (for example normal) with variance o"^ and zero mean. 

Let us describe the three terms in the right hand side. The first term is related to the 
marginal saving propensity of the agents, the second corresponds to the money transaction, 
and the last contains the effects of an open economy describing the market returns. Note 
that since debts are not allowed, the total amount of money in the system is increasing. 

This binary interaction model can be also related to some recent Lotka-Volterra type 
models jlj. A similar trade rule where r] and r/* have the same positive constant value 
has been considered recently in [^. In a closed economical system it is assumed that the 
total amount of money is conserved {r],rj^: = 0). This conservation law is reminiscent of 
analogous conservations which take place in kinetic theory. In such a situation equations 

correspond to a granular gas like model (or to a traffic flow model) where the stationary 
state is a Dirac delta centered in the average wealth (usually referred to as synchronized 
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traffic state in traffic flow modelling). Thus all agents will end up in the market with 
exactly the same amount of money. 

Let f{w^t) denote the distribution of money w E iR+ at time t > 0. By standard 
methods of kinetic theory, it is easy to show that the time evolution of / is driven by the 
following integro-differential equation of Boltzmann type, 

df 



dt 



(/?( 

J f{'w)f{'w:t:) - (3(w,w,)-^{w',wi)fi'w)f{w^)dW:t:dr]dr]^ (2) 



where Cw,'w<t:) are the pre-trading money that generates the couple {w,w^) after the 
interaction. In ((JJ J is jacobian of the transformation of {w,w^) into {w',w'^) and the 
kernel /? is related to the details of the binary interaction. 
We shall restrict here to a transition rate of the form 

P{w,w,)^{w',wi) = Kv)K^*)^{w' > 0)'i>{wl > 0), 

where ^{A) is the indicator function of the set A, and ^(•) is a symmetric probability 
density with zero mean and variance o"^. The rate function P[w',wi)^{w,wt) characterizes 
the effects of the open economy through the distribution of the random variables rj and rj^ 
and takes into account the hypothesis that no-debts are allowed. The above equation can 
be included in a more general settings where the trade rule has a more complex structure 
including, for example, risk, taxes and subsidies as described in 

We remark that, for general probability density //(•), the rate function /3 depends on 
the wealth variables {w,Wi:) through the indicator functions ^. This is analogous to what 
happens in the classical Boltzmann equation 0, where the rate function depends on the 
relative velocity. A simplified situation occurs when the random variables take values on 
the set (—(1 — 7), 1 — 7). In this case in fact, both w' > and > 0, and the kernel 
P does not depend on the wealth variables {w,w<t:). In this case the kinetic equation @ 
is the corresponding of the classical Boltzmann equation for Maxwell molecules i2j, which 
presents several mathematical simplifications. In all cases however, methods borrowed 
from kinetic theory of rarefied gas can be used to study the evolution of the function /. 



3 Main properties of the kinetic equation 

We will start our analysis by introducing some notations and by discussing the main 
properties of the kinetic equation 0. Let Q{f,f) denote the interaction integral. 



Q{f,f)(w) = 




(3) 

Let A^o the space of all probability measures in and by 

Mp = i^fieMo: j < +oo,p > o| , (4) 

the space of all Borel probability measures of finite momentum of order p, equipped with 
the topology of the weak convergence of the measures. 
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Let Ts{JR+)-, be the class of all real functions on iR_|_ such that ^(0) = ^'(0) = 0, and 
g{i^) ig Holder continuous of order 5, 

||.(-)|U = .upl''""W-<""<'">U^, (5) 

the integer m and the number < (5 < 1 are such that m + 5 = and (7*^™) denotes the 
m-th derivative of g. 

In the rest of the paper we will assume that the symmetric probability density ^(r/) 
which characterizes the transition rate belongs to A^2+a) for some a > 0. Moreover, to 
simplify the computations, we assume that this density is obtained from a given random 
variable Y with zero mean and unit variance, that belongs to M2+a- Thus, ^ of variance 
cr^ is the density of aY . By this assumption, we can easily obtain the dependence on a of 
the moments of /i. In fact, for any p > 2 such that the p-ih. moment of Y exists, 



/ \7]Y'n{r))dr) = E{\aY\P) = aPE{\Y\■ 
JlR 



By a weak solution of the initial value problem for equation (j^J, corresponding to 
the initial probability density fo{w) € Mp,p > 1 we shall mean any probability density 
/ G C^(iR+, A^p) satisfying the weak form of the equation 

^ [ ct>{w)f{w,t)dw = [Q{fJ)A) = 

"I J]R+ 

i IP(w,w,)-*{w',wi)f{'w)fiw^){(j){w')-(j){w))dw^dwdr]dr]^, (6) 
for t > and all (p £ J^p{lR+), and such that for all G J^p{]Rj^) 



lim / (f){w)f{w,t)dw= / (l){w)fo{w)dw. (7) 



We can alternatively use the symmetric form 

^ / f{w)4>{w)dw = \ I I /3(^,u-.)^KX)/(^)/(^* 



dt Jo 2 Jj[i2 Jjp(2 



+ 



(8) 

{4>{w') + (l){w'^) — (j){'w) — (j){'w^))d'w^ dwdrjdrj^. 

Existence of a weak solution to the initial value problem for equation Q can be eas- 
ily obtained by using methods first applied to the Boltzmann equation The main 
difference here is that the classical Boltzmann i/-theorem, which prevents formation of 
concentration, does not hold. 

^From © (or equivalently from ((H))) conservation of the total number of agents is 
obtained for (t){w) = 1 which represents the only conservation property satisfied by the 
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system. The choice (/'(w) = w is of particular interest since it gives the time evolution of 
the total wealth. We have 

^[ wf{w,t)dw= [ I [5(^^yj^)^(^,^^r)f{w)f\w^){'~i{w^-w) + riw)dw^dwdrjdrj^ = 

vKv)Kr]*) (l - ^{w' > 0)'^{w'^ > 0)) wf{w)f{w^)dw^dwdrjdrj^. (9) 

Since r/ < on the set where ^{w' > 0) = 0, equation shows that the total wealth 
is increasing, unless w' > and w'^ > for all exchanges (in such case the total wealth 
would be conserved). An upper bound for the total growth can be derived considering 
that 1 — 7 + ?7 > implies w' > (respectively 1 — 7 + r;* > implies w'^ > 0). Thus, 

1 - ^'(w' > 0)^'« > 0) < ^(1 - 7 + < 0) + ^-(1 - 7 + 7/* < 0). (10) 

Substituting this inequality into Q we obtain 

^/ wf{w,t)dw<-[ r//i(r?);u(r?*) (^-(1 - 7 + 7? < 0) 
dt Jm+ Jib? 

+ '^{l-'y + r]^ <0))dr]dT]^ wf{w)dw. (11) 

J1R+ 

By Markov inequality. 



/x(r/,)dr/,^'(l - 7 + r?* < 0) < _ 



while 

|??|/i(r?)(ir?*(l-7+r? < 0) < _^^y+o. l??P+°^(??)c^??*(l-7 + ?7 < 0) ^ 
Finally, 



2+a 



7) 



1+a ■ 



We remark that, if the kernel /3 does not depend on the wealth variables, the kinetic model 
preserves the total wealth. 

Similar bounds for moments of order higher than 2 can be obtained in a similar way. 
To this purpose, let us fix (j){w) = for some p > 2. Using the same trick as in Q we 
obtain 



^[ wPf{w,t)dw= [ I P(ui,w,)^(w',w',)f{w)f{w^){\w'\P -wP)dw^dwdr]dr]^ 
"I JlR+ Jm? J]R\ 

Kv)Kv*)f{'w)f{'^*){W\^ — w^)dw^,dwdr] d'q^ + 
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/ / fi{v)n{v*) (1 - *K > 0)^'K > 0)) {\w'\P - wP)f{w)f{w^)dw^dwdridri^ < 

/ / fi{r])iJ,{ri*)f{w)f{w:f){\w'\P — wP)dw*dwdri dri*+ 
Jib? J Ml 

[ [ fx{ri)fi{r]^)(l-^{w' >0)^{wi>0))wPf{w)f{w^)dw^dwdridri^. (13) 
JiH? J ml 



+ 

Using again Markov inequality we obtain 



/ / (1 - *K > 0)*K > 0)) wPf{w)f{w^)dw^dwdridri^ < 

Jm^ J Ml 

/ n{r])n{ri:,) - 7 + 77 < 0) + ^'(l - 7 + 77* < 0)) dr/ drj* / wPf{w)dw < 

JR^ Jm+ 



2 

Moreover, we can write 



\w'\P = wP +pwP ^ {w' - w) + ^p{p - l)\w\P '^{w'-wf, 
where, for some < 9 <1 

u, = 0u)' + 9)w. 



Hence, 



/ / M.)M*)/(»)/(».)(i»r-»')<i».<i»^,^,.= 




IJ'{ri)K'n*)f{'w)f{'^*){pw^ ^ {w' — w) + -p{p — l)\w\P '^{w' — w)'^)dw^dwdri drj^ < 
IM^ Jm% ^ 

\p{p-l) I I ii{r])ii{i^^)f{w)f{w^)\w\P~'^{w'-wfdw^dwdr]dr]^. (15) 

In fact, 

/ / M.)/.(*)/(w)/(w.)w'-'(w'-w)<iw.«,rf,.= 

Jm^ J Ml 

/ / ^{ri)ii{r]*)f{w)f{w*)vif~^{'~f{w* — w)+r]w)dwifdwdr}d'q* = 
Jm^ J Ml 

7 / f{w)f{w^)wP~^{w* — w)dw^dw < 0. (16) 
Jm^. 



The last bound follows by Holder inequality. In fact, 



f{w)f{wif)w'^ ^w^^dw^dw = I f{w)w^ ^dw / f{w*)w^dw^ < 
1R\ Jm+ J]R+ 

\ (P-1)/P / r \^^^ f 

f{w)w^dw I I / f{w^)w'^dw^, I = / f{w)w'^dw. 

1R+ J \J]R+ J J1R+ 



Finally, since 



\w 



we obtain 

fj.{ri)fi{r]^,)f{w)f{w^,)\w\^~'^{w' — w)'^dw^dwdri dr]^ < 
■ ['j'^{w^, — + Tffvj^ + 2r]w{vjt^ — w)\ dwtfdwdr] d-q^ < 



where 



CpApia,j) / wPf{w)dw, (17) 



^p(a,7) =^2(1 + 7^-2)+ + 272 + 27P + 27V^'-2. (18) 
Grouping all estimates together, we finally obtain the bound for the moments 



d 
di 



wPfiw, t) dw < Qp(p - l)cpApia, 7) + 2^^^^^) 1^ w^f{w)dw. (19) 



4 The continuous trading limit 

The previous analysis shows that in general it is quite difficult both to study in details the 
evolution of the wealth function, and to describe the asymptotic behavior. As is usual in 
kinetic theory, particular asymptotics of the equation result in simplified models (generally 
of Fokker-Planck type), for which it is relatively easier to find steady states, and to prove 
their stability. In order to give a physical basis on these asymptotics, it is relevant to 
discuss in some detail the interaction rule (^. To skip inessential difficulties, that as we 
will see later on do not change the structure of the limit equation, we suppose here that 
the random variables take values on the set (—(1 — 7), 1 — 7). As remarked in the previous 
section, in this case both tu' > and w'^, > 0, and the kernel [3 does not depend on the 
wealth variables {w,w.t). Let us denote by E{X) the mathematical expectation of the 
random variable X. Then the following properties follow from ^ 

E[w' + w'^] = w + w^, E[w' -w'^] = {l-2-i){w -w^). (20) 
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The first equality in H2U|) describes the property of mean conservation of wealth. The 
second refers to the tendency of the trade to decrease (in mean) the distance between 
wealths after the interaction. This tendency is a universal consequence of the rule in 
that it holds whatever distribution one assigns to ^, namely to the random variable which 
accounts for the effects of the market returns in an open economy. This universality is 
false for the first equality, which in general has to be substituted by the inequality 

E[w' + w'^]>w + w^. (21) 

Hence, in general the effects of the market returns in an open economy account for an 
increasing of the microscopic mean wealth. 

The second property in (|20() is analogous to the similar one that holds in a collision 
between molecules in a granular gas. There the quantity e = 27 is called "coefficient of 
restitution", and describes the peculiar fact that energy is dissipated j^H]- If we want to 
consider the situation in which most of the trades corresponds to a very small exchange 
of money (7 0), and at the same time we want to maintain both properties (|2()|) at a 
macroscopic level, we have to pretend that 

/ {w + w^)f{w)f{w^)dwdw^ = 2 wf{w)dw = 2m{t) (22) 
Jul. Jm+ 

remains constant, while 

/ {w-w^ff{w)f{w^)dwdw^ = Af{t) (23) 

is varying with time, and decays to zero when the market returns are not present (i.e. 
a = 0). 

When the kernel P does not depend on the wealth variables, © implies that m{t) = rriQ. 
Moreover, explicit computations show that Af{t) varies with law 



dAf{t) 
dt 



7(1-7) - ^ 



^/(t) + 2o-^m^ (24) 



Hence, when a = 0, Af{t) decays exponentially to zero. It is now evident that (|22j) is 
satisfied for any value of 7 and a, while H23() looses its meaning as 7 an cr tend to zero. Of 
course, one can rewrite (pi)) as 



dt 

Hence, if we set 



(l-7-^)^/(t) 



— m 

27 



(25) 



r = 7i, giw,T) = f{w,t), (26) 
which implies fo{w) = goiw), we obtain 
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Letting now both 7 — > and a ^ in such a way that = A, ()27() becomes in the 

hmit 

dAoir 



—^ = -{A-2X)Ag{r)+2Xm\ (28) 

This formal argument shows that the value of the ratio is of paramount importance 
to get asymptotics which maintain memory of the microscopic interactions. In particular, 
while for A < 2 ^g(r) converges to the finite value Ag = Am^/(2 — A) , Ag{T) diverges to 
infinity as time goes to infinity when A > 2. 

In the remaining of this section, we shall present a rigorous derivation of a Fokker- 
Planck model from the Boltzmann equation for the wealth function g{w,T), when both 
7^0 and a ^ in such a way that cT^/7 — > A. This derivation, which is similar to the 
quasi-elastic limit of granular gases, is of major relevance for the study of the asymptotic 
equilibrium states of the kinetic model. First, we show how the Fokker-Planck equation 
comes out for the simpler case of a kernel which does not depend on the wealth variables. 
Second, we extend the result to a general kernel. 

The scaled density g{v,T) = f{v,t) satisfies the weak form 

[ g4>dw = — [ [ ^{r])ji{r]^:)g{w)g{w^:){(]}{w') — (l){w))dw^.dwdrj drj^. (29) 
dr Jo 7 Jm^ Jul 

Given < 5 < a, let us set (p E J^2+s{lR+)- 
By ©, 

w' — w = 7(w* — w) + rjw. 
Then, if we use a second order Taylor expansion of (p around w 

4>{w') — (t){w) = {"fiw^ — w) + riw)(t)'{w) + ^(7(tf* — w) + rjw)'^ cj)" {w) , 
where, for some < < 1 

w = ew' + 9)w. 
Inserting this expansion in the collision operator, we get 

-1- I g(pdw = -[ f fi{7])fi{7]^)[{-f{w^ - w) +r]w)(l)'{w) + 
dr Jo 7 Jm^ J Ml 

+-(7(w* — w) + rjw)'^(l)'\w)]g{w^)g{w)dw^ dw drj drj^: + R{'j,a), (30) 



where 



i?(7,cj) = I [ fi{v)Kv*){l{w*-w)+Tjwy' 
" iM^ Jul 



27 

• [ct)"{w) - ct)"{w)) g{w.;,)g{w)dw^ dwdr]dr]^, 
Since (p £ !F2+&{JR+)-, and \w — w\ = 6\w' — w\ 

U"{w) - (P"{w)\ < U"\\s\w - wf < W'Wslw' - wf. (31) 
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Hence 



l«(7,^)l < 



27 



+ 



•|7(w* — w) + riw\^^^ g{w<^)g{w)dw^ dwd'qdr]^ 

By virtue of the inequality 

I7K -w)+ r]wf+^ < 4^+^ {\7w^\^+^ + |7u;|2+^ + \vw\^^^^ 
we finally obtain the bound 



\R{-f,a)\ <2 



1+2<5|| ,//| 



27^+^ + - / |r/|^+V('?)f^^ ) / -iy^+'sHdu- (32) 



1 



|2+(5, 



7 JiR 



K4 



Since is a probability density with zero mean and A7 variance, and fi belongs to M.2+aj 
for a > 6, 



[ \7]\'^-^^ fi{rj)dr] = E 

Jm 



X-fY 



2+5 



and E y\Y\^~^^j is bounded. Using this equality into (jHSJ one shows that R{'y, cr) converges 
to zero as 7 — > 0, if J^^^ w'^^^g{w, r) remains bounded at any fixed time r > 0, provided 
the same bound holds at time r = 0. By virtue of bound ([T^. taking = A7, we obtain 

d f 

— / nP''^^ g{w, t) dw < 
dr JM+ 

(A7)i+"/2 ' 



1/1 



w^g{w, T)dw. 



(33) 



iR4 



Since the lower order term in A2+5{\f^ , 7) is A7(l — 7)^, it follows that the boundedness 
of the moment holds independently of the value of 7. Therefore, at any fixed time r 

|i?(7,cT)|(T) < 2'+^'U"\\s [2^'^' + X^^'^"'E{\Yr')) Crj^ w^^'g,{w).dw (34) 

Hence, the remainder i?(7, cr) converges to zero as both 7 and a converge to zero, in such 
a way that cr^ = A7. Within the same scaling. 



lim — 

7^0 7 Jj[f Jjj2_ 



Kv)fJ'iil*)[ili'w* -w)+ r]w)(j)'{w) + 



+-(7(w* — w) + rjw) <j)' {w)]g{w^)g{w)dw^ dwdr/drj^. 



{m — w)(j)'{w) + —w (f)" {w) 



g{w)dw 



(35) 
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The right-hand side of (|35|) is nothing but the weak form of the Fokker-Planck equation 

(36) 



[uP'g) + -^{{w - m)g) . 



dr 2 dw'^ ^ ' dw 
The hmit Fokker-Planck equation can be rewritten as 

dg d , A , A 9 , 

a7 = 5^[(^-"^)+ 2^)5+2^9^(^5)]. 



(37) 



The general case of a wealth- variables depending rate function can be easily obtained 
from the previous computations. With respect to formula (|30j) we have to consider two 
more terms in the remainder. The first one comes out from the possibility to have negative 
wealth variables as outcome of the interaction. This term reads 



1 Jm' J Ml 



■ (l - -^{w' > 0)'^{wl > 0)) {4>{w') - (t>{w))g{w)g{w^)dw^dwdr] dr]^. (38) 

The integrand in 1)38(1 coincides with the integrand in the right-hand side of (|29|) , multiplied 
by the factor 1 — "^{w' > 0)^{w'^ > 0). Using ([T?H) shows that we have to integrate the 
probability density functions on a bounded domain, which gives a better behavior as 
7,(7 —> 0. In more details, while 



/ 



rf'^{rf) drj < 



1 



r]'^^"fi{v) drj < a^+^EiY 



2+a ■ 



'{|'?l>l-7} (1 -7)° y{|,,|>l-7} 

Thus, we can use the same expansion following formula ()29j) to conclude that Ri{'y,a) 
converges to zero as 7 ^ 0. 

The second remainder we have to take into account comes out from the fact that in 
general the total wealth m{t) is not constant in time. Thus we have the additional term 



^2(7,0-) 



1 



7 JiR2 \Jm+ 
On the other hand, thanks to (12 



w^g{w^:,T) dw^ — m ) g{w, T)(p' {w)dw^:dwdri dr]^. 

(39) 



/ 



w^g{w^,T) dw^, - m 



< 



exp 



7), 



(1-7) 



2+a 



m, 



which implies the convergence to zero of the remainder as 7 — > 0. Hence we proved 

Theorem 4.1 Let the probability density fo £ -Mp, where p = 2 + 5 for some 5 > 0, 
and let the symmetric random variable Y which characterizes the kernel have a density 
in M.2+a, with a > 6. Then, as 7 — > 0, cr — > m such a way that = A7 the weak 
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solution to the Boltzmann equation for the scaled density g^{v,T) = f{v,t), with t = ^t 
converges, up to extraction of a subsequence, to a probability density g{w^T). This density 
is a weak solution of the Fokker-Planck equation i':>b]) . and it is such that the mean wealth 
is conserved. 

The stationary state of the Fokker-Planck equation can be directly computed and, by 
assuming for simplicity 



it can be written as 



where 



m = f{vu,t)dw = l, 
Jr+ 



goo{w) = (40) 



r(/i) w^+t" 



Therefore the stationary distribution exhibits a Pareto power law tail for large w's. 

Note that this equation is essentially the same Fokker-Planck equation derived from a 
Lotka-Volterra interaction in j4, 2^ 2M ■ 



5 Numerical results 

In this section, we shall present some numerical test and we shall compare the stationary 
results obtained by using Monte Carlo simulation of the kinetic model with the stationary 
state of the Fokker-Planck model. We start from a situation where all individuals shares 
the same money. At each iteration, we choose randomly two individual. Then, the param- 
eter r] is choosen accordingly to a normal centered law (one for each individual) and the 
trade is performed if it is admissible i.e. when the new money of each individual remain 
positive. We report the results for the density f{w,t) = m(t) f {m{t)w , t)) where m{t) is 
the (variable in time) total amount of money in the system. Note that the density / is 
normalized to have the first moment equal to one. 

We use N = 2000 individuals and perform several iterations until a stationary state is 
reached. The distribution is then averaged over the next 250 iterations. 

We plot the distribution function for different values of (7, a^) = {(0.1, 0.2), (0.01, 0.02)}. 

The values for 7 and cr^ are such that we have a fixed A = 2 corresponding to a coeffi- 
cient jJL = 2 which is a realistic value for the income distribution of capitalistic economies 
^ 1^ • As prescribed from our theoretical analysis we observe that the equilibrium dis- 
tribution converges toward the function f^ as both 7 and o"^ go to 0, with (T^/7 = 2 (see 
Figures ^ and [2)) . 

The corresponding behaviors of the total amount of money obtained with N = 10000 
agents are given in Figure El Note that, as expected, asymptotically the exponential 
growth vanishes and the model preserves the total amount of money. 
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Figure 1: Asymptotic behavior of the Fokker-Planck model and the Boltzmann model for 
fi = 2.0, 7 = 0.1 and a = 0.2. Figure on the right is in loglog-scale. 

6 Conclusions 

We introduced and discussed a nonlinear kinetic model for a simple market economy, 
which is based on binary exchanges of money and speculative trading. We showed that 
at suitably large times, in presence of a large number of trades in which agents exchange 
a small amount of money, the nonlinear kinetic equation is well-approximated by a linear 
Fokker-Planck type equation, which admits a stationary steady state with Pareto tails. 
Convergence towards a similar steady state is shown numerically for the solution of the 
kinetic model after a suitable normalization which guarantees the conservation of the mean 
wealth. Our analysis enlightens both the range of validity of the Fokker-Planck equation 
H36|) . and why in the continuous trading limit the mean wealth, which is increasing for 
the kinetic model, remains constant. The present derivation takes advantages of the deep 
similarities between the trade rule and a molecular dissipative collision. However, it is 
expected that the asymptotic analysis discussed in this paper is not restricted to the trade 
rule (^, so that different microscopic interactions could be treated as well. It is clear that 
the analogy between the trade rule and a one-dimensional molecular dissipative collision 
suggests in a natural way the continuous trading asymptotic which is well-understood in 
kinetic theory with the name of quasi-elastic asymptotic procedure. Furthermore, the 
formation of overpopulated energy tails for large times in the kinetic model is in accord 
with the analogous result valid for the Boltzmann equation for a dissipative granular 
Maxwellian gas ^ EI- Pushing further these analogies may help to clarify many 
aspects of large-time behavior of market economies. 
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Figure 2: Asymptotic behavior of the Fokkcr-Planck model and the Boltzmami model for 
fi = 2.0, 7 = 0.01 and a = 0.02. Figure on the right is in loglog-scale. 
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